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Abstract 

In this paper we investigate asymptotic isotropization. We derive the asymptotic dynamics of 
spatially inhomogeneous cosmological models with a perfect fluid matter source and a positive cos- 
mological constant near the de Sitter equilibrium state at late times, and near the flat FL equilibrium 
state at early times. Our results show that there exists an open set of solutions approaching the de 
Sitter state at late times, consistent with the cosmic no-hair conjecture. On the other hand, solutions 
that approach the flat FL state at early times are special and admit a so-called isotropic initial singu- 
\Q • larity. For both classes of models the asymptotic expansion of the line element contains an arbitrary 

spatial metric at leading order, indicating asymptotic spatial inhomogeneity. We show, however, that 
C**) , in the asymptotic regimes this spatial inhomogeneity is significant only at super-horizon scales. 
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1 Introduction 

• i-H , 

In a recent paper we presented a general framework for analyzing the dynamics of generic spatially in- 
homogeneous cosmological models, referred to briefly as Go cosmologies (see Ref. [32], which we will 
henceforth refer to as Paper I). We employed Hubble-normalized scale- invariant variables which were 
defined within the orthonormal frame formalism, and led to the formulation of Einstein's field equations 
with a perfect fluid matter source as an autonomous system of evolution equations and constraints. In 
this paper we show that our framework can be used to derive asymptotic expansions for spatially inho- 
mogeneous cosmologies that undergo asymptotic isotropization at late times or at the initial singularity. 

It is part of the folklore of relativistic cosmology that all ever-expanding cosmological models with 
a positive cosmological constant asymptotically approach the de Sitter solution. A statement to this 
effect can be found in Gibbons and Hawking [10], p. 2739, in connection with a program to extend 
to cosmological event horizons certain ideas concerning thermodynamics and particle creation that had 
been applied to black hole event horizons. This conjecture, the so-called cosmic no-hair conjecture, 
was also made by Hawking and Moss [15], p. 36, after the introduction of the paradigm of cosmic 
inflation. Shortly thereafter, Starobinskii [30] took the first steps in describing the asymptotic structure 
of the metric for spatially inhomogeneous cosmologies with a perfect fluid matter source and an effective 
cosmological constant that approach the de Sitter solution at late times, showing that the line element 
contained 8 (the maximum possible number) free functions of the spatial coordinates. He commented 
that the asymptotic solution described "exponentially rapid local isotropization of the universe" , and 
coined the phrase "the cosmological constant is the best 'isotropizer'" . In the same year, Wald [33] 
gave a proof of the conjecture for ever-expanding spatially homogeneous (SH) cosmologies. Subsequently 
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Jensen and Stein-Schabes [17] extended Wald's proof to spatially inhomogeneous cosmologies, subject to 
the restriction that the 3-Ricci curvature scalar is negative. The significance of this result is not clear, 
however, for the following reason. The complete set of SH cosmologies can be divided into one invariant 
subset for which the 3-Ricci curvature scalar is never positive (Bianchi Type-I to Type- VIII) for which 
Wald's theorem applies, and another invariant subset containing Bianchi Type-IX and Kantowski-Sachs 
models, for which the dynamical behavior of the 3-Ricci curvature scalar is more complicated (some 
models approach the de Sitter solution while others re-collapse). Such a natural division into invariant 
subsets does not exist in the Go case. Hence, in Go cosmology the situation is still unclear as regards 
possible "cosmic no-hair theorems" . A variant of such a theorem, however, was proved very recently by 
Tchapnda and Rendall [31] for cosmologies with collisionless matter and positive cosmological constant 
that are either plane or hyperbolic-plane symmetrical. 

The motivation for much of the above work was provided by the idea of cosmic inflation, in that the 
de Sitter solution is regarded as describing the Universe in the late stages of an inflationary epoch. On 
the other hand, the de Sitter solution can also be regarded as a simple prototype of a cosmological model 
undergoing accelerated expansion during the present epoch. Indeed, recent observations have focused 
attention on such models. In particular, five years ago two independently working observational groups 
reported that the analysis of the type la supernova redshift data they had gathered suggested that in the 
present epoch the Universe may be in a state of accelerated expansion (see, e.g., Perlmutter et al [23], 
Schmidt et al [28], and Schmidt [27]). The acceleration could be due to either a veritable, repulsive 
cosmological constant, or some form of "dark energy"; see Carroll [6] and Ellis [9] for further details and 
references. 

In this paper, motivated by the above discussion, we consider Go cosmologies that are future asymp- 
totic to the de Sitter solution. We give a formal definition of this notion within the Hubble-normalized 
state space and then derive asymptotic expansions, valid as t — ► oo, for the Hubble-normalized variables. 

The second type of asymptotic isotropization, namely isotropization at the initial singularity, arises in 
connection with the notion of quiescent cosmology (cf., e.g., Barrow [2]), which provides an alternative to 
cosmic inflation. The idea is that, due to entropy considerations on a cosmological scale, a suitable initial 
condition for the Universe is that the Weyl curvature should be zero (or at least dynamically unimportant) 
at the initial singularity; this is the Weyl curvature hypothesis by Penrose [22], p. 630. This hypothesis 
leads to the notion of an isotropic initial singularity. Cosmological initial singularities of this type first 
arose as a special case in the general analysis of initial singularities performed by Lifshitz and Khalatnikov 
(LK hereafter) [19], p. 203. This type of initial singularity was also encountered in the work of Eardley, 
Liang and Sachs [8], p. 101. Subsequently, Goode and Wainwright [13] gave a formal definition of an 
isotropic initial singularity, using a conformally related metric, and derived various properties. Motivated 
by these ideas, we consider Go cosmologies that are past asymptotic to the flat FL solution and give a 
formal definition of this notion within the Hubble-normalized state space and then derive asymptotic 
expansions, valid as t — > — oo, for the Hubble-normalized variables. 

In this paper, we also extend the formalism of Paper I to cover geodesic null congruences and structures 
in connection with the issue of the formation of particle or event horizons (cf. the classic paper by 
Rindler [26]), thus naturally generalizing the work by Nilsson et al [21] from an SH context to a general 
spatially inhomogeneous setting. 

The plan of this paper is as follows. In Sec. 2, we present the Hubble-normalized evolution equations 
and constraints for Go cosmologies that arise from Einstein's field equations and the matter equations. In 
Sec. 3, we give the detailed asymptotic form of the Hubble-normalized variables for Go cosmologies that 
approach the de Sitter solution at late times, and discuss various features of this class of cosmological 
models. In Sec. 4, we present analogous results for Go cosmologies that approach the flat FL solution 
at early times. We conclude in Sec. 5 with a discussion of the analogies between the two classes of 
cosmological models that undergo asymptotic isotropization, and raise some issues for future study. The 
proofs of the results are given in the appendix. 

2 Evolution equations and constraints 

We consider spatially inhomogeneous cosmological models with a positive cosmological constant, A, and 
a perfect fluid matter source with a linear barotropic equation of state. We thus have 



= (7- 1)A . 



(1) 
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where fi is the total energy density of the fluid (assumed to be non-negative) and p its isotropic pressure, 
in the rest 3-spaces associated with the fluid 4-velocity vector field u, while 7 is a constant parameter. 
The range 

1 < 7 < 2 (2) 

is of particular physical interest, since it ensures that the perfect fluid satisfies the dominant and strong 
energy conditions and the causality requirement that the speed of sound should be less than that of 
light. The values 7 = 1 and 7=3 correspond to pressure-free matter ("dust") and incoherent radiation, 
respectively. 

We express the fluid 4-velocity vector field by 

«:=r(eo+«), (3) 

where e is a vorticity-free timelike reference congruence, v the peculiar velocity vector field of the fluid 
relative to the rest 3-spaces of e , and 

T := - 7f L= , v 2 := v a v a , (4) 

defines the usual Lorentz factor. To obtain an orthonormal frame, { e a } a =o.i,2,3, we supplement the 
vector field eo with an orthonormal spatial frame { e Q } Q= i j2 ,3 in the rest 3-spaces of e . The frame 
metric is then given by i] a b = diag [—1,1,1,1]. To convert the dynamical equations of the orthonormal 
frame formalism to a system of partial differential equations, it is necessary to introduce a set of local 
coordinates {x' 1 } M= o,i 1 2,3 = {*>#*}*=i,2,3- 



2.1 Hubble-normalized Einstein-Euler equations in separable volume gauge 

In Paper I, we employed the orthonormal frame formalism and introduced Hubble-normalized scale- 
invariant frame, connection and matter variables, as follows. 1 



{^o, da } 

{U a , A a , N aP , £ Q/3 , R a } 



{e ,e a }/H (5) 
{u a , a Q , n aP , a a p, fl a }/H (6) 



{SI, n A } := {/x, A}/(3ff 2 ) . (7) 
Expressing the Hubble-normalized frame vector fields with respect to a local coordinate basis leads to: 

d = AT 1 (d t - N* di) , d a = E a l , (8) 

where N and N l are viewed as coordinate gauge source functions, while the E a l are dependent variables 
that we refer to as the frame variables. 

The introduction of Hubble-normalized scale-invariant variables leaves H as the only dependent vari- 
able carrying a physical dimension, namely [ length The evolution equation for H decouples and is 
given by 

d Q H = -(q+l)H , (9) 

where q is the deceleration parameter, familiar from observational cosmology. The expression for q, which 
depends on the choice of temporal gauge, is given below. It is also necessary to introduce the spatial 
Hubble gradient r a , defined by 

r a :=-^d a H. (10) 

Throughout this paper, we will employ the separable volume gauge introduced in Paper I, which is 
characterized by the following choice of coordinate and frame gauge source functions: 

N l = 0, M=l => U a =r a ; (11) 



1 We use units such that Newton's gravitational constant G and the speed of light in vacuum c are given by 8ttG/c 2 = 1 
and c = 1. Then frame and connection variables have physical dimension [length] -1 , while curvature variables have physical 
dimension [ length ]~ 2 . 
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the choice N = 1 corresponds to setting the dimensional lapse function equal to the momentary Hubble 
radius, N = if -1 . With this choice of temporal gauge the deceleration parameter is given by 

q = 2Y? + ±G- + 1 [(? >1 -2) + (2- 1 )v 2 ]Cl-n A -\(d a -2A a )r a , (12) 

where S 2 = g (£ Q( 3£ a/3 ), and the scalars G± are defined by (see Eq. (6) in Paper I) 

G ± = l±( 1 -l)v 2 . (13) 

It should be noted that the temporal gauge choice (11) does not determine a unique family of spacelike 
3-surfaces S:{t = constant} and associated normal timelike reference congruence e$. In Sec. 4 [see 
Eq. (125) ] , we will use this remaining freedom to simplify the asymptotic expansions that we derive. There 
is also freedom in the choice of the spatial frame vectors and spatial coordinates. For the applications in 
this paper it is desirable to choose a Fermi-propagated spatial frame, i.e., one which satisfies 

R a = . (14) 

The remaining spatial gauge freedom is a time-independent rotation of the spatial frame vectors, 

~e a =Oj{x i )e p . (15) 

There is also the following freedom in the choice of spatial coordinates: 

x* (16) 

In the separable volume gauge with Fermi-propagated spatial frame, the Hubble-normalized state 
vector for Go cosmologies is given by 

X = (E a \r a ,A a ,N af3 ,Il af3 ,n,v a ,n A ) T . (17) 

The evolution equations and constraints that the components of X have to satisfy are given below. 2 

Evolution equations: 

d t Ej = (qSj'-Zj + e^R^E^ (18) 
d t r a = {q5 a p - 1 Lj + e aei p W)rp+a a q (19) 
d t A a = (q5 a p-^ a p + e a ^pR^)A^ + ±dp(Z af} + e a ^R->) (20) 

d t Z af3 = (q-2)Y, at3 ~2N^ a 1 N (3 ^ + N^ N^ -5'< ( - a (d 1 -r 1 )A r ^ 

+ e 75<Q [ (fl 7 - 2A y ) N^s + 2R 7 ] + (<F< a ()■ + A {a ) + 3 fiv< a i^ (21) 

G+ 

d t N af3 = (q5^ a s + 2^ a s + 2e'> s ^ R. ( ) N^ s -d 1 (e^ a ^\- 6^ a R f3) + 5 a<3 i? 7 ) (22) 

d t Q = --^v a d a n + G+ 1 [2G+q-{3 1 -2)-(2- 1 )v 2 - 7 (E^t/V 3 ) 
G+ 

- j(d a -2A a )v a +jv a d a \nG + ]fl (23) 

d t v a = - v 13 dpv a + 5 af3 dp In G+ - ^ - ^ (1 - v 2 ) S af3 (dp In Q - 2r ) 

7 

+ GZ 1 [ ( 7 - 1) (1 - v 2 ) (dpv' 3 ) - (2 - 7 ) v dp In G+ 

+ ll^Il ( 2 - 7 ) (1 - v 2 ) v S3 (dp In n - 2rp) + (3 7 - 4) (1 - v 2 ) 

7 

+ (2 - 7 ) (E^yV) + G - ( r P vf3 ) + i G + - 2 (7 - 1)] (Apvf 3 ) ] v a 
- Z a p v 13 + t a p 1 R v 7 - r a - v 2 A a + e af3 ~< Np S w 7 v 5 (24) 
d t fl A = 2 (q + 1) A . (25) 



2 These equations are obtained by imposing the restriction (11) on Eqs. (33)-(35), (38), (44), (47), (144) and (145) in 
Paper I. We keep the R a in the equations for future reference. 
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where q, G± and d a are given by Eqs. (12), (13) and (8). 
Constraints: 



= (C com ) i a p:=2(d [a -r [a -A [a )Ep ] i -e a psN^E J i (26) 

= (C G ) := 1 - Q k - S 2 - Q - Q A (27) 

= (C c ) a := dpZ af} + (2S a l3 - T, a p) r f3 - 3Ap £ Q/3 - e Q ^ Np S S 7 5 + 3 -2- Qv a (28) 

G+ 

= (Cj) a := {dp-rp){N af} + e a ^ A 1 )-2ApN af} (29) 

= (C A ) a := (0a-2r a )Q A . (30) 

Auxiliary 3-curvature variables: 

n k := -\(2d a -2r a -iA a )A a + \(N a pN a ' } )-^{N a a f (31) 
S Q/3 := - I e 75 (a (»| 7 | " r, 7 | - 2A W ) iV^ + ± (0 <Q - r <a ) Ag) + f 7V (c J N Ph - ± TV^ iV (a/3) (32) 

C a /3 := e 7<5 < Q (fl| 7 | - 2r\ 7 \ - A^\) Sp) S - 3iV( Q 7 5 (3 ) 7 + \ N^' S af} , (33) 



where S a p and fife are the tracefree part and trace part of the Hubble-normalized 3-Ricci curvature 
of a spacelike 3-surface S:{t = constant}, respectively, while C a p is the Hubble-normalized 3-Cotton- 
York tensor which provides information on the conformal curvature properties of a spacelike 3-surfacc 
S:{t = constant}. 

The 3-Ricci curvature variables S a p and satisfy the Hubble-normalized twice-contracted 3-Bianchi 
identity, given by 

= (dp - 2r p - 3Ap) S afj - e a ^ N ps S 7 S + | S af3 {dp - 2r fj ) Q k . (34) 



2.2 Energy-normalized null geodesies augmentation 

A few years ago, Nilsson et al [21] extended the Hubble-normalized equations for SH cosmology by adding 
the equations for an energy-normalized geodesic (timclikc or null) congruence. This strategy turned out 
to be very useful, both from an analytical and a numerical point of view. Wc therefore generalize that 
approach to the present Go cosmology setting, and this will subsequently allow us to investigate the 
connection between asymptotic silence of the gravitational field dynamics in relativistic cosmology and 
the formation of particle or event horizons (see Paper I, Sec. 4.1). 

Instead of being interested in the null geodesies emanating from a single given event, we will be 
interested in all possible null geodesies. We will thus not consider the problem of integrating a system of 
ordinary differential equations (associated with the former situation), but will instead consider all possible 
null geodesic vector fields (which, at a later stage, can be integrated to yield all possible associated flows). 
We therefore start by considering a vector field k a which is tangent to a geodesic null congruence. This 
vector field is thus governed by the equations 

k b W b k a = (35) 
k a k a = . (36) 

We now make a (3+l)-split of Eq. (35) with respect to a general orthonormal frame. Using the notation 
of Paper I and Ref. [21], this yields 

£e {£) = -k a (e a +u a ){£)-{H5 a [ ) + a a p)k a k f) (37) 
£e (k a ) = -k p ep{k a )-£{H5 a p + a a p)k fi -{kpk l3 )a a + {k fi ap)k a -e a p S n^k i:s k 1 

- £ 2 ii a + £ e a p 1 Q, 13 W , (38) 

while a (3+l)-split of Eq. (36) yields 

-£ 2 + {k a k a ) =0 . (39) 

Here £ := k° denotes the energy (frequency) of the particles that are traveling along the geodesic null 
congruence. 
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We now employ the Hubble-normalized gravitational field variables given in Eqs. (5) and (6) , and in 
addition introduce energy-normalized null vector variables as in Ref. [21] according to 

w 

K a := — . (40) 

With Eq. (39) we thus find 

K a K a = 1 , (41) 

where the K a correspond to the direction cosines of the null geodesies. 

Energy-normalization implies that the dimensional equation (37) for £ decouples and can be written 

as 

d £ = -(s + l)£, (42) 

where s is given by 

s := (E a/3 K a K^) - (t a - U a ) K a , (43) 

which is obtained from Eqs. (37), (40), (41) and (6). It is also necessary to introduce the spatial energy 
gradient t ai defined by 

t a :=-^d a £, (44) 

and governed by 

dot a = (qSj -Hj + e ai PRi)t + {d a -r a + U a ){s + l) (45) 
= {C t ) a := [ ^ {dp -rp-A p )-N^]t 7 . (46) 

The latter equations are obtained by choosing / = £ in the commutator equations (31) and (32) of Paper I, 
and making use of Eqs. (42) and (44) given above. Equations (45) and (46) constitute integrability 
conditions for Eqs. (42) and (44). 

We now write the evolution equation (38) in energy/Hubble-normalized form. 

8 K a = -K fi (d f3 -tf i -Ap)K a + { S 5 a f3 -Y, a p)Kf 3 -A a -e a p 5 N 1 s K f3 K-< 

- U a + e a (il R (i K'< . (47) 

Note that the contraction of Eq. (47) with K a vanishes identically on account of Eqs. (41) and (43). 
Let us now consider these equations in the separable volume gauge, specified by Eqs. (11). Then 

d t K a = -K f3 (dp-tp-A p )K a + (s5 a p-Z a p)K f} -A a -e a psN 1 s K l3 K~< 

-r a + e a Pl R fi if 7 (48) 

9tt a = (qSj-Xj + Ca/R^tp+daS, (49) 

where 

s = (Z aP K a KP) - (t a - r a ) K a . (50) 

Once the energy-normalized null vector variables K a have been obtained, one can solve for the associated 
geodesic null congruence by integrating the relation 

^-=E a ^K a 7 (51) 

subject to appropriate initial conditions. 

In this paper, we will use the above equations to determine asymptotic geodesic null structures and 
derive properties concerning particle and event horizons which are associated with asymptotic silence. 



3 de Sitter-like future asymptotics 

In this section we give the detailed asymptotic form of the Hubble-normalized variables for Go cosmologies 
that approach the de Sitter solution at late times. 
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3.1 G cosmologies future asymptotic to the de Sitter solution 

The de Sitter solution of Einstein's field equations describes a cosmological model with zero matter 
energy density and positive cosmological constant, which is undergoing exponential expansion. Employing 
a (3+l)-decomposition with intrinsically flat spacelike 3-surfaces, the line element and the associated 
timelike reference congruence are given by 

ds 2 =-dT 2 + £% c 2 VW (dx 2 + dy 2 + dz 2 ) , (52) 
e = cV , (53) 

where A > is the cosmological constant, and T is clock time. It is convenient to choose the unit of 
[length] as follows: 

to = yfl ■ (54) 

The Hubble scalar is constant and positive, and is given by 

H =\f^ = i o 1 - (55) 
It follows from Eq. (9) that the deceleration parameter is constant and negative, being given by 

q = - 1 . (56) 
We also find it convenient to use a conformal time coordinate r), defined by 

n := e~ T / e ° . (57) 

The line element is then cast into the form 

ds 2 = f 2 r,- 2 (- dr, 2 + dx 2 + dy 2 + dz 2 ) . (58) 

Relative to the natural orthonormal frame associated with the line element (52), the Hubble- normalized 
variables are 

E a l = e"* 5 a l , r a = , (59) 

S Q/3 = 0, A a =0, N aP = 0, (60) 

o = o, fi A = 1 , ( 61 ) 

where the dynamical time coordinate t is introduced via 

t=y, (62) 

or, equivalently, 

tj = e"* . (63) 

We note that the volume density V is given by 

V = i% e 3t , (64) 

so that the separable volume gauge conditions are satisfied. 

Observe that the frame variables satisfy lim^oo E a l = 0. Thus, the de Sitter solution is described by 
an orbit in the Hubble-normalized state space that is future asymptotic to an equilibrium point on the 
silent boundary, E a l — 0. 

Motivated by Eqs. (59)-(61), we say that a Go cosmology is future asymptotic to the de Sitter solution 
if the following limits are satisfied: 

lim {E a \r a ) T = (65) 

t — >oo 

lim(£ Q/3 ,yr,7V Q/3 ) T = (66) 

t— too 

lim Q = , lim fl A = 1 . (67) 

t^oo t^-oo 
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We note that for the de Sitter solution the peculiar velocity variables v a are unrestricted, since the 
Codacci constraint (28) (which may be viewed as determining the v a algebraically provided that f! ^ 0) 
yields no information. We will show that if the limits (65)-(67) are satisfied, together with certain 
technical restrictions, then the evolution equations determine the asymptotic form of the v a as t — > oo. 
Wald's theorem [33], mentioned in the introduction, asserts that any SH cosmology with a positive 
cosmological constant, except those of Kantowski-Sachs and Bianchi Type-IX which may re-collapse, is 
future asymptotic to the de Sitter solution, in the above sense. As discussed in the introduction, at present 
it is not known whether Wald's theorem can be generalized to Go cosmologies. We now consider the class 
of Go cosmologies which are future asymptotic to the de Sitter solution, and present the asymptotic form 
of the Hubble- normalized variables as t — > oo. 



3.2 Asymptotic expansions 

In App. B, we prove that the de Sitter solution is asymptotically stable, and derive the asymptotic form of 
the Hubble-normalized variables for Go cosmologies that are future asymptotic to the de Sitter solution, 
in the sense of Eqs. (65)-(67). It is necessary to impose certain restrictions on the spatial derivatives; 
these are given in App. B. 

We now give the asymptotic expansions, using the convention that "hatted" coefficients depend only 
on the spatial coordinates: 



(E a \A a ,N a P) J 



(£ Q i ,i Q ,7V a/3 ) e-* + 0(e- 3t ) 

- \ (E a l d^k) c- 3t + 0(c-(W + e" 5t ) 

- 3S al3 c- 2t + t af} e~ 3t + 0(e' 4t ) 
l-Cl k c- 2t + 0{e~ 3 ^ +e' 4t ) 

f2e- 4t +0(c- 5t ) 



for 1 < 7 < | 
for 7 = | 



f2e- 4t +0(e- 5t +e' 



A n (37-4) i t . 



for | < 7 < 2 



= < 



$a e (3 7 -4)t 



+ 0(e" 3 *) 
0(e~ 



.3(37-4)* 



1 — v = exp 



+ 0(e~*) 



for 7 = 1 
for 1 < 7 < 
for 7 = | 



(1 - v 2 ) + 0(e-* + e" 



(2-7) 

The coefficients A a and N al3 are determined by E a l according to 

A a = E p l E [ J d 3 E^ 



(2-7 



(68) 

(69) 
(70) 
(71) 

(72) 



(73) 



for | < 7 < 2 .(74) 



where the matrix E a i is the inverse of E a % : 

E a t = S a t 



E a i Ej = Si j . 



(75) 
(76) 

(77) 



The coefficients &k and S a p in Eqs. (69)-(71) are the leading order coefficients in the asymptotic expan- 
sions of the 3-Ricci curvature variables, 



n k = c - 2t Ci k + o{ c - 2t ) 

S afi = e~ 2t [S a(i + 0{e~ 2t ) 

They are determined by E a l , A a and N a p according to: 3 

O fc = -\{2E a l d t -?,A a )A a + \{N a pN a P)~±{N a a ) 2 



(78) 
(79) 

(80) 

S af} = - | e<\ a {E h \% - 2A h{ ) N p)s + | E (a % A (j) + f N p)l - \ N (a[j) . (81) 



3 These equations are the "hatted" versions of Eqs. (31) and (32), with r a set to zero. 
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The constraint (Cc) a at order e 4 * provides the restriction for v a : 

= (Ej d l - iA p ) t af3 - e a ^ N ps £y 5 + C a + Q a , (82) 

where 

r , r - • „ a ( 3 7^« a for 1< 7 < |. 

C = i Q = i 12Qv a /{3 + vpv' 3 ) for 7=3. (83) 

^ 7 ' (3&v a for|< 7 <2. 

The constraints (Ca) q and (Cq) have been used in App. B to give the coefficients in Eqs. (69) and (71), 
respectively. 4 



3.3 Features of asymptotic to de Sitter G cosmologies 

In this section we investigate certain features of asymptotic to dc Sitter Go cosmologies, as described by 
the asymptotic expansions (68)-(74). Our first goal is to show that these expansions represent a general 
class of perfect fluid Go cosmologies. We accomplish this by showing that the expansions (68)-(74) 
contain 8 freely specifiable functions of the spatial coordinates x % , the same number that appear in the 
initial data for a general cosmological solution of Einstein's field equations with a perfect fluid matter 
source (see LK, p. 188). 



Essential arbitrary functions 

The nine coefficients E a l in the expansions (68)-(74) can be chosen as suitably differentiable arbitrary 
functions of the x l , which then successively determine A a , N a p, Clk and S a p via Eqs. (75), (76), (80) 
and (81). The five shear rate coefficients S Q( 3 and the coefficient VL in the density parameter can also be 
chosen arbitrarily, while the peculiar velocity coefficients v a are determined algebraically by the Codacci 
constraint (82), employing Eq. (83). As mentioned in Sec. 2, the separable volume gauge (11) leaves a 
freedom in the choice of the 1-parameter family of spacelike 3-surfaces S:{t = constant}. A coordinate 
transformation of the form 5 

t = t + <p(x i ) + 0(e- 2t ) , 5* = ^-!^ g ij djipe- 2t + 0(e- 3t ) , (84) 

where g 1 ^ = £q 2 S a/3 E a l , preserves the separable volume gauge to leading order. This freedom can be 
used to fix Q, although we do not choose to do so in general. Three of the E a l can be eliminated by the 
frame rotations (15), and three by the spatial coordinate transformations (16), leaving three essentially 
arbitrary functions in the E a l . The asymptotic expansions (68)-(74) thus contain 8 essentially arbitrary 
functions of the spatial coordinates. That is, they represent a general class of perfect fluid G cosmologies 
that are future asymptotic to the de Sitter solution. 



Metric expansion 

We now derive an asymptotic expansion for the spacetime metric in order to relate our results to the 
work of other researchers, in particular that of Starobinskii [30]. We find that the metric expansion sheds 
light on the nature of large-scale spatial inhomogeneity in Go cosmologies that are future asymptotic to 
the de Sitter solution. 

In the separable volume gauge, the line element has the form 

ds 2 = - R- 2 dt 2 + gij dx l dx j , (85) 

where [see Paper I, Eq. (163)] 

gij = H- 2 5 a0 E a t E p s . (86) 

4 Note that up to (and including) order e — 2t , the present expansion also satisfies the gauge conditions defining the 
constant mean curvature (r„ = 0) and synchronous (U a = 0) temporal gauges, respectively. 

5 See Starobinskii [30] for the analogous synchronous gauge-preserving transformation. The function ip(x l ) can be inter- 
preted in terms of a boost of the original frame by the boost function w a {t, x l ) via w a = — <5 Q ' 3 (Ep 1 dup) c~' + 0(c~ 2t ). 



3 DE SITTER-LIKE FUTURE ASYMPTOTICS 



10 



Equations (68)-(74) are substituted into Eq. (12) to give an improved expansion for q, and then into 



Eq. (18) to give an expansion for d t E a l , which can be integrated. In the case 7 < |, we obtain 



73 



(87) 



The components E a i of the corresponding Hubblc-normalizcd 1-forms form the inverse of the matrix 
E a l , and are given by 



E a i = e* 



s a f3 + U3S a f3 + n k s c 



-it 



-3t 



Sl6 a p e- arft + C(e- 4t ) 



with E^j and E^' related by Eqs. (77). It also follows from Eq. (9) that the Hubble scalar H has the 
expansion 

H = £- 1 \l + ±Q k c- 2t + ±£10-^ + <D(c- 4t )] . (89) 
We then substitute Eqs. (88) and (89) into Eq. (86) to obtain the following expansion for the 3-metric gtj : 

0(e- 4t )" , 



9ij 



5 a a E a i E ,3 1 



gij + 3Sij c 2t - 

S a n E a i & 



2 y.. P -3t 

3 ^3 e 



E a/3 E c 



(90) 
(91) 



where 

Qij &af3 E i E^ j , (S^j O a fj j-j , . — , — , ( .l . , 

Observe that the free functions E a l determine the coefficients g\j in the leading order term in the expan- 
sion through Eqs. (77) and (91), and also determine the second term §ij, while the free functions 
enter at a higher order. Note that up to order C(e _4t ), the matter coefficient O only enters the spacetime 
metric (85) through H in Eq. (89). It can be shown that the expansions (89) and (90) are also valid for 
7 satisfying | < 7 < 2, with the difference that the matter coefficient enters into the 0(e _4t ) term in 
Eq. (89). 

Our results provide a confirmation of the ansatz for the metric expansion given by Starobinskil [30] , 
Eq. (2). Some of the details differ, however, due to the fact that Starobinskil employs the synchronous 
gauge. We have established consistency with his results by performing a coordinate transformation of 
the form 

i = T = j H- 1 dt + . . . , £' = &' + .... (92) 

We have also derived Starobinskii's metric expansion directly in the synchronous gauge, using our inte- 
gration method. 

The most striking feature of the expansion (90) is the fact that the leading order coefficient in the 
expansion is an arbitrary 3-metric, g^, whereas the spacelikc 3-surfaces S:{t = constant} in the de 
Sitter line element (52) are intrinsically flat. The generality of the leading order coefficient g^ in the 
metric expansion raises an apparent paradox — a Go cosmology that is future asymptotic to the spatially 
homogeneous and isotropic de Sitter solution in the sense of the definition in Subsec. 3.1 can exhibit 
substantial spatial inhomogeneity. We will show that this paradox is resolved by the existence of event 
horizons. 



Event horizons 

In cosmology an event horizon for a fundamental observer can be thought of as constituted by the 
observer's past light cone at t = 00 (see Hawking and Ellis [14], p. 129). In order to establish the 
existence of event horizons, we need to determine the asymptotic form of geodesic null congruences. 
These are governed by Eq. (51), with the K a determined by Eqs. (48)~(50). 
In App. D we show that any null geodesic has the asymptotic form 

- a4 - K a E a i (x J 00 ) e~* + 0( C - 2t ) , (93) 

for constants x 1 ^ and K a . For given x l x , the family of null geodesies then form an event horizon for 
the fundamental observer whose worldline is x' 1 — x^. The spatial distance from the observer to her/his 
event horizon is given by 

*<" -jf \/^¥? d » • < 94 > 
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where y l — y l (s) (with < s < 1) describes a spacelike geodesic from to x l (t), for fixed t. In the 
asymptotic regime, the spacelike geodesic is basically a straight line, and is approximated by the null 
geodesic (93) projected onto a spacelike 3-surface S:{t = constant}. It follows from Eq. (93) that 

d H (t) = *o + 0(e-*), (95) 

where £q is given in Eq. (54). 

One can always introduce local coordinates at x l — x 1 ^ such that (see e.g., Schutz [29], p. 156) 

g tf (0 = £l Sa + 'f (x k - at) (*' -xl) + ... ■ (96) 

For points within the event horizon the approximation cjij{x m ) w l\ Sij becomes increasingly accurate as 
t — > oo, due to Eq. (93). In other words, within the event horizon of a particular fundamental observer, 
the spacetime metric asymptotically approaches the de Sitter metric. 6 Since gij is general, however, the 
above approximation cannot be done simultaneously at all points, reflecting the spatial inhomogeneity of 
the Go cosmology at super-horizon scales. We emphasize that the spatial inhomogeneity does not diminish 
as t — > oo — it is the observer who sees successively smaller portions of that spatial inhomogeneity. 



Radiation bifurcation 

A noteworthy feature of the asymptotic expansion (68)-(74) is that for the physically important case of 
incoherent radiation (7 = |), a bifurcation occurs that affects the peculiar velocity v of the perfect fluid 
relative to the spacelike 3-surfaces S:{t — constant}. 

If 1 < 7 < |, the components v a tend to zero, and as a result the fluid 4- velocity vector field u 
asymptotically coincides with the vorticity-free timelike reference congruence e^. For incoherent radiation, 
the v a do not tend to zero in general, and v 2 can approach any value between and 1. If | < 7 < 2, 
v 2 tends to the value 1 as t — > 00, which means that the speed of the fluid relative to the fundamental 
observers comoving with eg approaches the speed of light. 

It is of considerable interest to investigate the limit behavior of the fluid kinematical variables, as 
perceived by the fundamental observers. Using the boost transformation laws provided in App. F, we 
find that the future asymptotic limits for the Hubble-normalized fluid kinematical scalars (defined by 
Eqs. (270)-(272) in App. F) are 

Urn (U 2 , S 2 , W 2 ) = (0, 0, 0) for 1 < 7 < f . (97) 

lim {U 2 , S 2 , W 2 ) = \ i w 2 , 0, 1 for 7 = \. (98) 

lim (U 2 , S 2 , W 2 ) = [3( 7 -l) 2 , H 3 7-4) 2 , 0] for § < 7 < 2. (99) 

This extends previous results given by Goliath and Ellis [11] for some special SH cases with cosmological 
constant to the general Go case (see also Goliath and Nilsson [12] and Raychaudhuri and Modak [24]). The 
context of spatial inhomogeneity stresses the isotropization issue. We have to ask ourselves isotropization 
along which timelike congruence? In the present context there is no isotropization along the timclikc 
reference congruence eo when 7 > | , but this does not say anything about whether the models isotropize 
along the fluid 4-vclocity vector field in this case, which, perhaps, is the physically best motivated 
congruence to consider. However, to determine whether this is the case or not takes us outside the scope 
of the present paper, and so we leave this issue for future studies. We note that when 7 < |, the fluid 
4-velocity vector field asymptotically coincides with the timelike reference congruence and, hence, one 
does have isotropization in this case. 



4 Flat FL-like past asymptotics 

In this section we give the detailed asymptotic form of the Hubble-normalized variables for Go cosmologies 
that approach the flat FL solution at early times. 

6 Boucher and Gibbons [4] have clarified this matter in a different way, by transforming the de Sitter metric to a static 
form, valid within the event horizon of a particular fundamental observer. 
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4.1 G cosmologies past asymptotic to the flat FL solution 

The flat FL solution of Einstein's field equations describes a cosmological model with perfect fluid matter 
source, whose density parameter has the constant value fi = 1. The line element and the fluid 4- velocity 
vector field are given by 

/ rp \ 4/(37) 

d.s 2 = -dT 2 + ^ 2 ( — J (dx 2 + dy 2 + dz 2 ) , (100) 
u = d T , (101) 
employing clock time T as one of the local coordinates. The matter energy density and pressure are 

[l =±^T-\ p=( 7 -l) M . (102) 

The Hubble scalar is 

H=^T~ 1 , (103) 
37 

while the deceleration parameter is given by 

9=1(37-2); (104) 

the latter is thus constant and positive. We find it convenient to use a conformal time coordinate r], 
defined by 

, T y3 7 -2)/3 7 3 7 T 

' % = (*7=2) 4' (1 ° 5) 

with the unit of [length] £ and the constant T chosen to satisfy r]o = 1, i.e., 

37 



T . (106) 



(37-2) 

The line element is then cast into the form 

d S 2 = £ 2 ,,4/(37-2) (_ d?? 2 + dj ,2 + dy 2 + ^ (10?) 

and the Hubble scalar is 

H = Ato" 1 ^- 3 ^- 2 ) . (108) 
We introduce the dynamical time coordinate t via 

/ rp \ 2/(37) 

y -«'• (io9 » 

or, equivalcntly, 

^ = e J(37-2)t. (110) 

Relative to the natural orthonormal frame associated with the line element (100), the Hubble-normalizcd 
variables are [using Eq. (106)] 

Ej = \{^-2)e^-^5j (111) 
£ Q/3 = 0, A a =0, N a p = 0, (112) 

= 1, O\ = 0. (113) 

We note that the volume density V is given by 

V=^e 3 *, (114) 

so that the separable volume gauge conditions are satisfied. 



4 FLAT FL-LIKE PAST ASYMPTOTICS 



13 



Observe that the frame variables satisfy limt^-oo E a l = 0. Thus, the fiat FL solution is described 
by an orbit in the Hubble-normalized state space that is past asymptotic to an equilibrium point on the 
silent boundary, E a l = 0. 

Motivated by Eqs. (Ill)— (113), we say that a Go cosmology is past asymptotic to the flat FL solution 
if the following limits are satisfied: 



lim (E a \r a ) T = 

t—>- — oo 

lim {V a f,,A a ,N a f,) T = 



lim = 1, 



lim flf 



. 



(115) 
(116) 
(117) 



We note that for the flat FL solution the peculiar velocity variables v a are zero. We will show that 
if the limits (115)-(117) are satisfied, together with certain technical restrictions, then the constraints 
imply that the v a tend to zero as t — > — oo. We now consider the class of Go cosmologies which are past 
asymptotic to the flat FL solution, and present the asymptotic form of the Hubble-normalized variables 
as t — > — oo. 



4.2 Asymptotic expansions 

In App. C, we derive the asymptotic form of the Hubble-normalized variables for Go cosmologies that are 
past asymptotic to the flat FL solution, in the sense of Eqs. (115)— (117). It is necessary, as in Subsec. 3.2, 
to impose certain restrictions on the spatial derivatives; these are given in App. C. 

We now give the asymptotic expansions. For brevity we use n instead of e^ 37-2 )* [see Eq. (110)]. 



(E a i ,A a 1 N a > 3 y 



E n 



s a/3 = _ 



A a ,N af) Y r? + C(r? 3 ) 
VL k )rf + 0{rf + 
S al3 r 1 2 +0{r ] 4 ) 



\ (E a * diClk) r, 3 + 0( V 5 + if+WTT,) 
6 



(37 + 2) 



A 



V 



2+ 



(37-2) 



1-fW + W + rf + WT2j) 
' "a/3 ( £)_o_ "\ ^3 i fi^fnJy i „3+" 



■ 5 a " (V diSi k ) ^ + o( v b + n^Ts^s) 



(37 + 2) 

where H is the leading order coefficient in the asymptotic expansion of H: 



H = Hrf 



[i+ 



(37-2)1 



l + ±n k n 2 +0(r, 4 



2+7 



") 



(118) 
(119) 
(120) 

(121) 

(122) 
(123) 

(124) 



In the derivation in App. C, we used the freedom to re-define the 1-parameter family of spacelike 3-surfaces 
S:{t — constant} while preserving the separable volume gauge to set 



H 



37 



1 o ' 



(125) 



so that H = constant > 0. In addition, Eq. (10) and the constraints (Cq) and (Cc) a were used in App. C 
to give the coefficients in Eqs. (119), (122) and (123), respectively. 7 

As in Subsec. 3.2 before, the coefficients A a , N a p, Cl k and S a p are determined by the E a % according to 
Eqs. (75), (76), (80) and (81). Likewise, the 3-Ricci curvature variables have the asymptotic expansions 



(lk+0(r, 2 ) 



(126) 
(127) 



7 Note that up to (and including) order rj 2 , the present expansion also satisfies the gauge conditions defining the constant 
mean curvature (r a = 0), synchronous (U a = 0) and fluid-comoving (v a = 0) temporal gauges, respectively. 
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4.3 Features of asymptotic to fiat FL G cosmologies 

In this section we investigate certain features of the class of Go cosmologies that are past asymptotic to 
the flat FL solution. One of our goals is to show that this class of cosmologies admit an isotropic initial 
singularity in the sense of Goode and Wainwright [13] (see also Newman [20] and Anguige and Tod [1]). 



Essential arbitrary functions 

Of the position-dependent coefficients in Eqs. (118)-(123), only the nine E a l can be chosen arbitrarily, 
as suitably differcntiable functions of the x l . However, three of them can be eliminated by the frame 
rotations (15), and three by the spatial coordinate transformations (16). Thus, the asymptotic expansion 
contains 3 essentially arbitrary functions of the spatial coordinates, as compared with 8 for a general class 
of perfect fluid Go cosmologies. These models thus form a set of measure zero in the Hubble-normalized 
state space. 



Metric expansion 



We now derive an asymptotic expansion for the spacetime metric in order to relate our results to the 
work of other researchers, in particular that of LK and Goode and Wainwright [13]. We find that the 
metric expansion sheds light on the nature of large-scale spatial inhomogeneity in Go cosmologies that 
are past asymptotic to the flat FL solution. 

In the separable volume gauge, the line element has the form given in Eq. (85). Equation (18) gives 
an expansion for d t E a l , which can be integrated to give 



^(3 7 -2)(3 7 + 2) 

Following the steps taken in Subsec. 3.3, we obtain the metric expansion 



m = V 



4/(37-2) 



12 



(3 7 - 2)(3 7 + 2) 



Sn v 2 + o( v 4 + n 



2+ 



(37-2) 



where 



9ij = H 2 <W E a i E 



3 > 



Sij — H 2 S a Q E a i E 



(128) 



(129) 



(130) 



Our result provides a rigorous confirmation of the ansatz for the metric expansion made by LK, Eq. (4.1), 
for the case of incoherent radiation. Some of the details differ, however, due to the fact that LK employ 
the synchronous gauge. We have established consistency with their results by performing a coordinate 
transformation of the form given in Eqs. (92). We have also derived LK's metric expansion directly in 
the synchronous gauge, using our integration method. 



Isotropic initial singularities 

We briefly digress to introduce the notion of an isotropic initial singularity; see Ref. [13]. A cosmological 
model 8 is said to admit an isotropic initial singularity if the physical spacetime metric g is conformal to 
an unphysical metric g: 

9 = X 2 9, (131) 

where \ = x{v) i s a differentiable function of a time coordinate r\ and satisfies x(0) = 0, while the 
conformal metric g is regular on the spacelike 3-surface S:{n = 0}. The conformal factor x(v) 1S a ls° 
required to satisfy 

lim *M = oo (132) 

lim xWx^y 2 



3 We assume that the matter source is a perfect fluid satisfying Eqs. (1) and (2). 
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It follows from Eqs. (85), (124) and (129) that the spacetime metric for the class of Go cosmologies that 
are past asymptotic to the flat FL solution does satisfy Eqs. (132) and (133) with 

X = W /(37 " 2) , (134) 

where r] is the conformal time coordinate. 
Particle horizons 

The generality of the leading order coefficient gij in the metric expansion (129), which should be compared 
with Eq. (90), raises an apparent paradox — a Go cosmology that is past asymptotic to the spatially 
homogeneous and isotropic flat FL solution in the sense of the definition in Subsec. 4.1 can exhibit 
substantial spatial inhomogencity. We will show that this paradox is resolved by the existence of particle 
horizons. 

In order to establish the existence of particle horizons, we need to determine the asymptotic form of 
geodesic null congruences. These are governed by Eq. (51), with the K a determined by Eqs. (48)-(50). 
In App. E we show that any past-oriented null geodesic has the following asymptotic form as rj — > 0: 

j(v) = *bb + J^Y } ^4Vbb) V + 0(r, 2 ) , (135) 

for constants x BB and K a . The constants x BB give the point of termination of the null geodesic at the 
singularity r\ = 0. Thus the past null cone at time rj for the fundamental observer whose worldline is 
x l = x will intersect the singularity r\ — 0, thereby defining a particle horizon, which is formed by the 
fundamental worldlines given by x l — x BB . The spatial distance from the observer to her/his particle 
horizon at time 77 is given by 

% ds < (i36) 

where y l = y l (s) (with < s < 1) describes a spacelike geodesic from x BB to x % (77), for fixed r\. As in the 
derivation of Eq. (95), it follows from Eq. (135) that 

1+ 



dn(v) =V 1+ T&=V [e + O(r])} , (137) 

where £0 is given in Eq. (106). 

One can introduce local coordinates at x l = x BB such that 

= ll S tj + (** ~ *bb) (*' ~ 4b) + • • ■ • (138) 

For points within the particle horizon the approximation gij{x m ) w £q Sij becomes increasingly accurate 
as r] — > 0, due to Eq. (135). In other words, within the particle horizon of a particular fundamental 
observer, the spacetime metric is close to the fiat FL metric. Since is general, however, the above 
approximation cannot be done simultaneously at all points, reflecting the spatial inhomogeneity of the 
Go cosmology at super-horizon scales. We emphasize that the spatial inhomogeneity does not diminish 
— it is the observer who sees successively smaller portions of that spatial inhomogeneity. 



5 Concluding remarks 

In this paper, we have given a unified analysis, within the framework of the Hubblc-normalizcd state 
space, of Go cosmologies that undergo asymptotic isotropization, either at late times (asymptotic to the 
de Sitter solution, with £ — > 00) or near the initial singularity (asymptotic to the flat FL solution, with 
£ -» 0). 

The analysis reveals a number of common features as regards the asymptotic behavior in the two 
classes, as well as a number of significant differences, which emerge most clearly if we write the asymptotic 
expansions in Subsec 3.2 in terms of the conformal time coordinate of the de Sitter line element, given 
by Eq. (57). We note that the length scale factor is given by 



£(r 1 )=£ aV - 1 , 



(139) 



5 CONCLUDING REMARKS 



16 



for the de Sitter solution, and by 

*(»7) = W /(37_2) , (140) 
for the flat FL solution. Thus, for both asymptotes rj — > + in the asymptotic regimes under consideration. 

Firstly, the shear rate S Q( g and anisotropic 3-Ricci curvature S a p have the same leading asymptotic 
dependence on rj as rj — > + : 

E a ^c4H 2 + 0(>) 3+5 ) (141) 
5^-5 Q/3 r? 2 + 0( ?7 3 + 5 ) , (142) 

where 5 > 0, and C = — 3 in the de Sitter case while C = — 6/(37 + 2) in the flat FL case. 
Secondly, there is a common asymptotic form for the line element, namely 

ds 2 =£ 2 (r ] )[-dT] 2 +£o 2 g ij dx i dx j +0{ V 2 )] , (143) 

as rj — > + , where £(77) is the length scale factor for the de Sitter solution or for the flat FL solution. 
Likewise, the Hubble scalar has the common form 

H = H asyiap ( V ) [i + ±fW + eV +5 )] , (144) 

where 5 > 0, and H asymp (r/) is the Hubble scalar for the de Sitter solution or for the flat FL solution, 
respectively [see Eqs. (55), (89), (108), (124) and (125)]. Here Cl^ is the leading order coefficient in the 
expansion for the 3-Ricci curvature scalar £l fe , which has a common form for the two classes, namely, 

fi fc =fW + 0fa 4 ). (145) 

The key feature of the asymptotic line element (143) is the presence of the arbitrary 3-metric gij(x k ). 
We have shown that the spatial inhomogencity that it generates is significant only at super-horizon scales 
(the event horizon for Go cosmologies that are future asymptotic to the de Sitter solution and the particle 
horizon for Go cosmologies that are past asymptotic to the flat FL solution). 

The common asymptotic features of the two classes can be understood to some extent within the 
context of the Hubble-normalized state space. For each class, the asymptotic state is described by an 
equilibrium point on the silent boundary, 

E a l = , 

[see Paper I, Eq. (75)], which is associated with an isotropic and conformally flat solution of Einstein's 
field equations. Use of a conformal time coordinate rj, which tends to zero on the silent boundary for 
both classes, leads to common asymptotic decay rates for its various Hubble-normalized quantities. The 
asymptotically silent dynamics (E a l — > 0) is also responsible for the existence of an event horizon and a 
particle horizon, as shown in Subsecs. 3.2 and 4.2. 

The differences between the two classes originate in the asymptotic behavior of the deceleration 
parameter q, which is negative in the de Sitter case (q — > —1) and positive in the flat FL case (q — > 
5(37 — 2)). This difference affects the higher-order terms in the asymptotic expansion of the shear 
rate £ Q/ g. In the de Sitter case, there is a term £ Q/ 3?7 3 , where the £ Q( 3 are freely specifiable functions, 
while in the flat FL case, this term does not appear. 9 The absence of the £ a| g-term, of course, accounts 
for the difference in the number of freely specifiable functions in the two cases. This difference is also 
reflected in the asymptotic expansions for the Hubble-normalized electric and magnetic parts of the Weyl 
curvature [see Paper I, Eqs. (152) and (153)]. For the electric Weyl curvature we have in the de Sitter 
case 

£ a p = It^rf + Otf) , (146) 

and in the flat FL case, 

£ a p = j^ Y) S a(3 r i 2 + 0{r ! i ). (147) 
For the magnetic Weyl curvature we have in both cases 

Hap = - C aP rf + Otf) , (148) 

9 The second shear rate mode in this case is a growing mode (into the past), and the assumption of the isotropization 
requires the corresponding coefficient to be zero. 
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where C a p is given by the "hatted" version of Eq. (33), with r a set to zero. 

The difference between the dc Sitter and the flat FL cases as regards the number of freely specifiable 
functions can be interpreted in a dynamical systems context. We have already noted that both the de 
Sitter solution (52) and the flat FL solution (100) determine equilibrium points on the silent boundary in 
the Hubble-normalized state space. It follows that the orbits of Go cosmologies that are future asymptotic 
to the de Sitter solution (respectively, past asymptotic to the flat FL solution) are asymptotic to these 
equilibrium points. The asymptotic expansion contains 8 arbitrary functions (the maximum possible 
number) in the de Sitter case, which confirms that the de Sitter equilibrium point is asymptotically 
stable, as proved in App. B (see Stage 0). On the other hand, the presence of only 3 arbitrary functions 
in the flat FL case implies that the flat FL equilibrium point has both a stable and an unstable manifold 
(into the past), and our asymptotic solutions describe the stable manifold. 

A major accomplishment of the present paper is to provide derivations of the detailed asymptotic 
properties of Go cosmologies that undergo asymptotic isotropization, instead of simply making ad hoc 
assumptions about the form of the metric, as has been done previously. Our derivations do, however, 
depend in a significant way on the assumption that the spatial partial derivatives of the Hubble-normalized 
variables remain bounded in the asymptotic regimes. It is this assumption that enables us to rely 
exclusively on analytic methods for systems of ODE. 

It would certainly be desirable to attempt to use methods from the theory of PDE to prove the 
validity of the assumption of bounded spatial derivatives for cosmological models that undergo asymptotic 
isotropization. Recent experience with G2 cosmologies, which we now describe, gives some cause for 
optimism. Firstly, it is known that in G2 cosmologies spatial partial derivatives can diverge on approach 
to the initial singularity, through the creation of so-called Gowdy spikes (see Paper I, Subsec. 4.1, and 
Ref. [25]). This spatial structure is created, however, by the local instability of the Kasner solutions. Since 
these solutions do not play a role in the present context (the approach to an isotropic initial singularity) , 
we do not expect spatial structure of this nature to develop. Secondly, at late times, support for our 
assumption is provided by numerical simulations 10 of G2 cosmologies, which do not show the development 
of large spatial derivatives. 

A second limitation of our analysis is that it is local in nature in that we restrict our considerations 
to some open subset of the integral curves of the timclike reference congruence e^. This restriction is 
inevitable, however, since in a Go cosmology all timelines do not necessarily share the same asymptotic 
evolution. For example, a Go cosmology with positive cosmological constant may approach the de Sitter 
solution along some timelines but may undergo collapse to a future singularity along others. 

Another matter that requires comment is the choice of temporal gauge. In the present paper we have 
chosen to work with the separable volume gauge, defined by Eqs. (11), because we have found it to be 
a convenient computational gauge. On the other hand, we have been able to transform our asymptotic 
expansions to the synchronous gauge, and find that the leading order terms, e.g., Eqs. (141)-(148), are 
unchanged. This calculation thus provides evidence for the gauge robustness of our results. 

In conclusion, it would be of interest to relate our results to perturbation analyses of the de Sitter 
solution (see for example Barrow [3] and Bruni et al [5]) and of the flat FL solution. We have not been 
able to clarify the relationship to our satisfaction, particularly in the case of the de Sitter solution, and 
so we leave this matter for future investigation. In addition, other possible applications of the Hubble- 
normalized state space framework for Go cosmologies presented here that suggest themselves naturally are 
establishing firm links with mainstream issues in observational cosmology such as the cosmic background 
radiation, peculiar motions of galaxies, or gravitational lensing. 
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A Asymptotic results for ordinary differential equations 

Proposition 1 

Given 

d t M a[i = AJMrf + B aP , (149) 
where M a p, A a p and B af} are Cartesian tensors. 11 If 

lim A af3 (t) = -k8 af3 , k > , (150) 

t — >oo 

and 

B aP = 0(e- lt ) , l>0, (151) 

then for any given e > 0, 

M af} = 0(e(- m+e ^) , (152) 

where m = min(fc, I). 
Proof: Introduce 

M = (M af} M al3 ) 1/2 , M a& = MS a/3 , where {S al} S al3 ) 1/2 = 1 . (153) 

It follows from Eqs. (149) that 

d t M = (S af3 A a ^S lP ) M + S af3 B a0 . (154) 
Now Eq. (150) implies that for all e > there exists a to such that 

\S af3 {A a i + kS a iS lP \<e , (155) 
for all t > t , and Eq. (151) implies there exists a C > such that 

\S af> B a0 \ < Ce~ lt . (156) 

It then follows that 

d t M < (-k + e)M + Ce~ u . (157) 
Multiplying by e' fc_e ^ and integrating from t = to to t, yields 



M < 



Ce -lto 

M - 



{-k + e + l) 



e (-fc+ e )(t-to) + _±1 e -«(*-*o) ; (i5 8 ) 

( K 6 I ) 



from which the result follows. □ 

Proposition 2 

Given 

d t M afj = AJM^p + B aP , (159) 
where M a p, A a p and B a p are Cartesian tensors. If 

lim A a e{t) = k5 afi , k > , (160) 

B a p = 0(e~ H ), l>0, (161) 

and M a p is bounded, then 

M aP = 0(e- lt ) . (162) 



11 We will also need a vector version of this result, i.e., the DE is of the form 

d t M a = A a /3 M' 3 + B a . 
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Proof: As in the proof of Proposition 1, we obtain Eqs. (154) and Eq. (155) with k replaced by —k, 
and Eq. (156). It follows from these equations that 

d t M > M(k - e) - Ce~ lt . (163) 

Multiply by q-^-*) 1 and integrate from t to oo, knowing that M is bounded, to obtain 

M ^(^7T7) e " Zt ' < 164 > 

and the result follows. □ 

Proposition 3 

Given 

d t M = - kM + B , (165) 
where M and B are m x n matrices and k > is constant. If B = 0(e~ lt ), I > k, then 

M = M c- kt + 0(e- lt ) . (166) 
Proof: Rewrite the differential equation (DE) as d t (e kt M) = e kt B, and integrate. □ 

B Derivation of asymptotically de Sitter future dynamics 

In this appendix we derive the asymptotic expansions (68)-(74) for Go cosmologies that arc future 
asymptotic to the de Sitter solution in the sense of Eqs. (65)-(67). For ease of reference, we list the 
complete set of restrictions below. 

Al: E a \ A a , N aP , r ai S Q/ 3, and J1a — 1 are sufficiently small at some initial time which we can take 
to be t = 0. 

A2: The partial derivatives of any dimensionless variable V are bounded as t — > oo. In addition, if V 
tends to zero, then diV = 0(\\V\\) as t — > oo, for fixed x 1 , where ||V|| is the magnitude of V. 12 

A3: If V = Ve~ nt + C(e~ mt ) as t -» oo, where m > n > 0, then d t V = diV e~ nt + 0{e~ mt ) as t -» oo. 
Stage 

We first establish asymptotic stability of the de Sitter solution. The evolution equations for the variables 
in Al arc of the form 

d t X = AX + f{t,X) , (167) 

where X is the vector containing the variables in Al, A is a constant diagonal real- valued matrix with 
negative eigenvalues, while f(t,X) contains products of components in X and d a X . Assumptions Al 
and A2 imply that f(t,X) is continuous for all t > and f{t,X) = o(\\X\\) as ||X|| 0. Then 
Theorem 1.1 in Coddington and Levinson [7], p. 314, implies that X = is asymptotically stable. Thus, 
we obtain the property 

PI: Ea, A a , N a p, r a , Y, a p, 0, Oa - 1 — » as t — » oo, for fixed x l . 

We organize the derivation of asymptotic expansions into two stages, in which we repeatedly use the 
evolution equations and the constraints (Cq) and (C\) a . In Stage 1, we obtain order estimates for all 
variables, primarily using Propositions 1 and 2 in App. A. In Stage 2, we use the results from Stage 1 in 
conjunction with Proposition 3 to obtain the explicit asymptotic expansions. The assumptions A2 and 
A3 are needed throughout in order to bound the spatial derivative terms in the evolution equations and 
in the constraints (Cq) and (CtC) a . 



2 For example, \\EJW = S a ^E a i Ep i , for each i, and ||S a£ 
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Stage 1 

The deceleration parameter q, given by Eq. (12), plays a dominant role in the evolution equations. Its 
limiting value as t — > oo follows from PI and A2: 

lim q = - 1 . (168) 

t — >oo 

We begin by considering E a l . Using PI and Eq. (168), the evolution equation (18) for E a l assumes the 
form of Eqs. (149) and (150) with k = 1 and B aj 3 — 0. It follows from Proposition 1 that for any given 
e E > 0, 

E a l = e>(e(- 1+e ^*) , (169) 

as t — > oo. 

We can now use Eqs. (168), (169) and PI, A2 to conclude that the evolution equation (20) for A a 
assumes the form of Eqs. (149) and (150) with k = 1 and i = 1 — eg. It follows from Proposition 1 that 
for any given €a > £e, 

A a = C(e ( ~ 1+£A)t ) , (170) 

as t — > oo. 

At this stage we have two different epsilons in Eqs. (169) and (170), with ea > £e- Equation (169) 
is equally valid, however, if we replace ce by ca, and in the interests of simplicity we choose to do this, 
dropping the subscripts on the epsilons. We will also make this e-simplification in subsequent steps. 

In a similar way, Eq. (22) leads to 13 

N a @ = 0(e<- 1+£ >*) . (171) 
Next, we derive intermediate results for O and r a . For any given S > 0, PI and A2 imply that 14 

— r'd.M • . (172) 

for t sufficiently large. Then, for any given e > 0, Eqs. (23) and (168)-(172) give 



^(3 + « 2 ) + £ 



d t n < 

for t sufficiently large. The inequalities 1 < 7 < 2 imply that 



n , (173) 



-7^(3 + « 2 )<-3. (174) 

It follows from Eqs. (173) and (174) by integrating that 

fl = 0(c ( - 3+ ^) . (175) 

The constraint (C\) a , in conjunction with PI, A2 and Eq. (169), implies 

r a = 0(e(- 1+£ )*) . (176) 

With Eqs. (175) and (176), we now have enough information to derive an asymptotic expression for S™' 3 . 
Indeed, Eq. (21) assumes the form of Eqs. (149) and (150) with k = 3 and I = 2 — e, which implies, using 
Proposition 1, that 

= C( e (- 2+e )*) . (177) 

It follows from Eqs. (169) (171), (176) and A2 that tt k = C(c^ 2+e )*), which in turn implies, using 
Eqs. (175), (177) and the constraint (Cg), that 

fi A - 1 = 0(e(- 2+e )*) . (178) 



13 Here and in the future, we will omit the qualifier t —* 00. 

14 Here we use A2 to write di^l = 0(Q), and also use the fact that E a l — > to write \E a l \ < e, for t sufficiently large. 
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By using A2 and the constraint (C\) a , we can now conclude that 

r a = Oi^- 3 ^ 1 ) . (179) 

It remains to derive the decay rates of f2 and v a . We now make use of the evolution equation for the 
scalar v = ^v a v a , which can be derived from Eq. (24). The equation for v 2 is in fact given in Paper I; 
see Eq. (146). Using the preceding asymptotic expressions and A2, the equation for v assumes the form 

8 t v = l-{l-v 2 ) (37- A)v + g , (180) 

where g = 0(e^~ 1+e " >t ). Equation (180) is in fact an asymptotically autonomous DE, a scalar version of 
the vector DE discussed by Horwood et al in Rcf. [16]. We can thus use Theorem B.l in Ref. [16], p. 15, 
to obtain the limiting behavior of v as t — > oo. The domain D for v is the interval < v < 1, and the 
conditions Hi and Hi in Ref. [16] are satisfied. It is straightforward to verify that the solutions of the 
autonomous DE 

d t v=-^-(l-v 2 )(3 1 -A)v , (181) 



with initial conditions in D, satisfy 



if 1 < 7 < 



4 



lim v = < - ' 3 ' (182) 

I 1 if | < 7 < 2. V ' 



t— *oo 



It then follows from Theorem B.l in Ref. [16] that any solution of Eq. (180) with initial condition in D 
satisfies Eq. (182). 

If 1 < 7 < §, Eq. (180) assumes the form of Eqs. (149) and (150) with k = -(37 - 4) and I = 1 - e. 
Proposition 1 then implies 

v = 0{e^- 4+e ^) . (183) 
If I < 7 < 2, we write the evolution equation for v [cf. Paper I, Eq. (146)] in the form 

d t {l-v 2 ) = -v a d a (l-v 2 )- -I- [(3 7 -4)i; 2 +0(e ( - 1+e) *)l (1 - v 2 ) . (184) 

PI and A2 imply that for any given 5 > 0, 15 

-v a d a (l -v 2 ) < 5(1- v 2 ) . (185) 

Then Eq. (184) assumes the form of Eqs. (149) and (150) with k = 2 and B afi = 0. Proposition 1 

implies 

1 - v 2 = (D(c [ - 2 7^W +e]t ) . (186) 

For 7 = |, Eq. (24) gives 

d t v a = 0(e^ 1+ ^) , (187) 

which implies 

v a = v a + 0(e ( ' 1+e)t ) . (188) 
Using Eq. (182), the inequality (173) can be strengthened to read 

t n<{^ + ^ forl< 7 <!, (18Q) 
\ (-4 + e) n for § < 7 < 2, V ' 

where the e has been redefined. If 7 = |, Eq. (173) simplifies directly to 

d t tt < (-4 + e) Q . (190) 



Inequalities (189) and (190) now give 



'e>( e (-3 7 +e)t) forl<7<§. 
0(e(- 4 + £ )*) for I < 7 < 2. 



3 Sincc 1 — v 2 — » in this case, we can use A2 to write — v 2 ) = 0(1 — v 2 ), and proceed as in footnote 14. 
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Stage 2 

In the second stage we use the asymptotic decay rates found in Stage 1 to successively write the evolution 
equations in a form to which Proposition 3 can be applied. 

We begin by considering E a l . The decay rates in Stage 1 and Eq. (12) gives 

q = - 1 + 0(e(- 2+e >*) . (192) 

The evolution equation (18) for E a l can now be written in the form 

d t E a l = -E a i +gJ , (193) 

where g a l — 0(e^ 3+c '*). This equation, in conjunction with Proposition 3, implies that 

E a l = E a l e~* + 0(e(- 3 + e )*) . (194) 

In a similar way the evolution equations (20) and (22) lead to 

A a = A a c- t + 0{e ( -~ 3+e '> t ) (195) 
N aP = N af) e-* + 0{c { - 3+e)t ) . (196) 

We now consider the evolution equation (23) for O. For 1 < 7 < |, we obtain 

d t n = - 3~fQ + g , (197) 

where g = 0(e( 37 ~ 8+c )*), with flv 2 being the dominant term in g. For | < 7 < 2, we obtain 

d t Q = - 40 + g , (198) 

where g = 0{e ( -- 5+£ '> t + e ( " 4 " 2l ^T +e )*), withd Q 0, nd a v a , (A a v a ) Q and Q {l-v 2 ) being the dominant 
terms in g. For 7 = |, we obtain 

d t n = -4O + 0(c(- 5+£ )*) . (199) 

Applying Proposition 3 yields 



'(] e -37f +C )( e (37-8+ e )t) forl< 7 <§. 



ft c- 4 * + 0(e(- 5 + e )* + e { - 4 - 2 l^W +t)t ) for § < 7 < 2. (200) 
fte~ 4t + 0(e^- 5 +^) for 7=|. 

At this stage, the evolution equation for E a/3 assumes the form 

d t T, af) = - 3E Q/3 - 3S af) c- 2t + g a!i , (201) 

where g a ^ = 0(c(~ 4+e )'). By imitating the proof of Proposition 3, we conclude that 

At this stage (Cq) and A3 give 

Oa - 1 = - ft k e- 2t + 0(e~ 3 ^ + e~ 4t ) , (203) 

while (C\) a and A3 lead to 

r a = -\ (Ej dift k ) e- 3t + C( e -( 1+3 '') t + c" 5 *) . (204) 
We now consider the evolution equation (24) for v a . For 7 = 1, we obtain 

d t v a = - v a + C(e(- 3+e )*) . (205) 

Proposition 3 implies 

v a = v a e- 1 + 0( e (- 3+e >*) . (206) 
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For 1 < 7 < |, we obtain 

d t v a = (37 - 4) v a + g a , (207) 

where g a = 0(e"* + e 3 ( 3 T- 4+e )*), with d Q lnfi and v 2 v a being the dominant terms in g a . Proposition 3 
implies 

v a = v a e ( 3 T" 4 )* + 0( e -* + e 3(37-4+«)*) . (208) 
For 7 = |, we already have the result: 

i;" = -0" + 0(e(- 1+e)t ) . (209) 
For | < 7 < 2, we obtain from Eq. (184) 

ft (1 - v 2 ) = 2 fci) (1 - t, 2 ) + ff , (210) 

where ff = 0(e (_1 " 2 ^T +e)t + e (_4 ^T +£)t ), with (l-v 2 )fi a w Q , (l-i; 2 )A a -u a and (l--y 2 ) 2 being the 
dominant terms in g. Proposition 3 implies 

1 - , 2 = (1 - v 2 ) e~ 2 ^ + O^- 1 - 2 ^^ + e(- 4 ^ + ^*) . (211) 

Using the results obtained in Stage 2, we can repeat Stage 2 to eliminate the epsilons in the O-terms 
as follows. Equation (192) now reads q = — 1 + 0(c~ 2t ), and the function g Q 4 in Eq. (193) satisfies 
g a l = 0(e -3 *); similarly for the other variables. 

The position-dependent coefficients E a l , A a and N a p are required to satisfy the constraints 

= 2 (E [a i dj - A [a ) Erf - e a0S N 5 ~< E^ (212) 
= En 1 d t (N af3 + e a ^ Ay) - 2A/3 N af3 , (213) 

which arise, respectively, from (C com ) 4 a/3 and (Cj) Q , at order e -2 *. Equation (212) can be solved to express 
the coefficients A a and N a p in terms of the coefficients E^ and their inverse E a i. These expressions, 
given by Eqs. (75) and (76), satisfy Eq. (213) identically. The constraint (Cc)" provides a restriction on 
v a , given by Eqs. (82) and (83). □ 



C Derivation of asymptotically flat-FL past dynamics 

In this appendix we derive the asymptotic expansions (118)-(123) for Go cosmologies that are past 
asymptotic to the flat FL solution in the sense of Eqs. (115)— (117). For ease of reference, we list the 
complete set of restrictions below. 

Al: Ej, A a , N afj , r a , S Q/3 , l-O, fi A -> as t -> -co, for fixed x\ 

A2: The partial derivatives of any dimensionless variable V are bounded as t — > — oo. In addition, if V 
tends to zero, then diV = 0{\\V\\) as t — > — oo, for fixed x\ where ||V|| is the magnitude of V. 

A3: If V = Ve nt + 0(e mt ) as t -» -oo, where m > n > 0, then d l V = d t V e nt + 0{e mt ) as t -» -oo. 

We organize the derivation into two stages, in which we repeatedly use the evolution equations, (10) 
and the constraints (Cq) and (Cc) Q - In Stage 1, we obtain order estimates for all variables, primarily 
using Propositions 1 and 2 in App. A. In Stage 2, we use the results in Stage 1 in conjunction with 
Proposition 3 to obtain the explicit asymptotic expansions. The assumptions A2 and A3 are needed 
throughout in order to bound the spatial derivative terms in the evolution equations, in (10), and in the 
constraints (Cq) and (Cc) a ■ 
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Stage 1 

First, Al, A2 and (C c ) a imply that {~//G+) Ctv a -> 0. But we have fi -> 1. Thus we have 

lim v a = . (214) 

t^ — oo 

It follows from Al, A2 and (12) that 

lim q = \ (37 - 2) . (215) 

We now introduce r = — t, for convenience in using Propositions 1-3. 

We begin by considering E a l . Using Al and Eq. (215), the evolution equation (18) for E a % assumes 
the form of Eqs. (149) and (150) with k = \ (37 — 2) and B a p = 0. It follows from Proposition 1 that 
for any given eg > 0, 

E a l = 0(e [ -^- 2)+€E]T ) , (216) 

as t — > 00. 

We can now use Eqs. (215), (169) and Al, A2 to conclude that the evolution equation (20) for A a 
assumes the form of Eqs. (149) and (150) with k = 5 (37 — 2) and / = | (37 — 2) — e^. It follows from 
Proposition 1 that for any given €a > e_E, 

A a = 0(e [ -^- 2)+eA]T ) , (217) 

as r — > 00. We make the e-simplification as in Subscc. B. 
In a similar way, Eqs. (22), (19) and (25) lead to 16 

N a P = 0{c [ -^~<- 2)+e]T ) (218) 

r a = (D(e [ -^- 2)+e]T ) (219) 
n A = 0(c^ + ^ T ) . (220) 

With Eqs. (216)-(219) and A2, the constraint (C c ) a gives 

v a = 0{e [ -^- 2)+i]T ) . (221) 

Equation (21) now assumes the form of Eqs. (159) and (160) with k = | (2 — 7) and I = (37 — 2) — e. 
Since S Q/ 3 is bounded asr^oo (it is assumed to tend to zero), Proposition 2 implies 

E a/3 = e>( c [-( 3 7-2)+ e ]r) _ ^ 222 ^ 

Lastly, the constraint (Cq) implies 

1 - n = 0(c[-( 3 t- 2 )+ £ ] t ) . (223) 

Stage 2 

In the second stage we use the asymptotic decay rates found in Stage 1 to successively write the evolution 
equations in a form to which Proposition 3 can be applied. 

We begin by considering E a l . The decay rates in Stage 1 and Eq. (12) give 

q = \ (3 7 - 2) + 0( C [-(3T-2)+ £ ]r ) _ (224) 

The evolution equation (18) for E a l can now be written in the form 

= - i (3 7 - 2) Ej + g a * , (225) 

16 Here and in the future, we will omit the qualifier r — » 00. 
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where g a l = C(c' 2( 37 2 )+ e l T ). This equation, in conjunction with Proposition 3, implies that 



E a l = E a l c -l^-<- 2 ) T + 0{c [ -^-<-V+^ T ) . 
In a similar way the evolution equations (9), (25), (20), (22) and (19) lead to 

H = i?C (37/2)T + 0( e [(37/2)-(3 7 -2)+e]T) 

A a = A a c^ { ^^ 2)T + 0(c h l (37_2)+£lT ) 

N afi = i ya/3 e -i(3 7 -2)r +c)( - e [-|(37-2)+ £ ]r^ 



where Cl A = A/(3# 2 ), and f a = - (E a l d l H)/H from (10). 
At this stage, (Cc) a and A3 give 



37 



■f a e" 2 



1 (3 7 -2)r + £) ^[-|(3 7 -2)+e]r- ) 



while (Cg) and A3 lead to 



0=1- Cl kC -^- 2 ^ + ( c [-2(3 7 -2)+e]r + & 



■3 7 t\ 



(226) 

(227) 
(228) 

(229) 

(230) 

(231) 



(232) 



(233) 



We now use the freedom to re-define the 1-parameter family of spacelike 3-surfaces S:{t = constant} 
while preserving the separable volume gauge 17 to set H = (2/37) T ~ , with T a positive real-valued 
constant. As a result, 

f a = , 



and v a = 0(c [ -2^->- 2 1+^). 

The evolution equation (21) for assumes the form 

d T T, af) = |(2-7) S a/3 + 35 Q/3 c - (37 - 2)r + g ap , 
where g afj = 0( c [-2(3 7 -2)+ e ]r^ We rewritc Eq (235) as 

,), ( ( ,-5( 2 -^S a ^ =35 a/3 e- (1+ i 7)r + ff a/3 e-| (2 - 7)7 



and integrate from r to 00 to obtain 

E a/3 = - ■ 



. 5«/3 e -(3 7 -2)r + 0( e [-2(3 7 -2)+e]r) _ 



(37 + 2) 

Lastly, we derive the leading order coefficients for r a and v a . First, 

g = § (3 7 - 2) - \ (3 7 - 2) il k C ^- 2 )- + 0(c [-2(3 7 -2) + e]r + c -3 7 r } 

Evolution equation (9) for H then yields 

H = H C^^ T [ 1 + \ Ofc e"^" 2 ^ + 0( c [-2(3 7 -2)+ e] r + e -3 7 r } " 

Eq. (10) and the constraint (Cc) a then successively give 

\ (Ej di(l k ) + 0( c [- 2 ^- 2 H £ l T + c- 3 ^)" 



r«=c-2( 3 ^ 



(234) 



(235) 



(236) 



(237) 



(238) 



(239) 



(240) 



17 A transformation of the local coordinates of the form t = t+tp(x % )+0(e^ 3 ~ f 2 *> t ),x l = x % + ^ 2 <f J 9jy? e^ 37 2 )'- 

0( e 2(3 7 -2)t _|_ e 3 7 t^ w j[[ p reserve the separable volume gauge. 
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and 

y a =e -f(37-2)T r fia + ( e [-2(37-2)+ e ]T + e -37Tjl > ( 2 41) 

where 

« Q = <5 Q/3 (S/s* W + ^ + (Ej d t S a(3 ~ 3A fj S a ? - e a ^ N ps £/) . (242) 

The twice-contracted 3-Bianchi identity (34) gives 

= Ep* diS° fi - SA P S afi - e a ^ N 0S + ± S afj {Ej ^O fe ) . (243) 

This simplifies v a to 

fi ° = (37T2) r ' 3( ^^ ) ' (244) 
Using the results obtained in Stage 2, we can repeat Stage 2 to eliminate the epsilons in the 0-terms, 
as in App. B. 

As in App. B, the position-dependent coefficients E a % , A a and N a p are required to satisfy the con- 
straints (212) and (213), which arise, respectively, from {C com ) l a f} and (Cj) Q , at order e~^ 37_2 ' r . Equa- 
tion (212) can be solved to express the coefficients A a and N a p in terms of the coefficients E a l and their 
inverse E a i. These expressions, given by Eqs. (75) and (76), satisfy Eq. (213) identically. □ 

D Asymptotic expansions for event horizons 

Equation (44) and expansion (68) imply that t a = C(e _t ) as t — > oo, subject to the assumption that £ 
satisfies A2. Equations (43) and (42) successively imply that 

s = 0(e-*) , £ = 0(e-*) . (245) 

With these results, the evolution equation (48) for K a yields d t K a = 0(e~ t ), which can be integrated 
to give 

K a = k a + 0{c- t ) . (246) 
With expansions (68) and (246), Eq. (51) is integrated to give 

x i (t)=x i O0 +0(e- t ) . (247) 
The coefficients E a l [x 3 ' (t)] are then Taylor-expanded: 

E a i [x\t)}=E a i {xi ) + 0{e- t ) . (248) 
With expansions (68), (246) and (248), Eq. (51) is integrated to give an improved expansion 

x\t) = xl - K a E a i (x> 00 )e- t + 0(e~ 2t ) . (249) 

E Asymptotic expansions for particle horizons 

Equation (44) and expansion (118) imply that t a = 0(rj) as r\ — > 0, subject to the assumption that l/£ 
satisfies A2. Equations (43) and (42) successively imply that 

s = 0{r]) , l/£ = 0{if /{ ^- 2) ) . (250) 

With these results the evolution equation (48) for K a yields 

K a = K a + 0{r}) . (251) 

With expansions (118) and (251), Eq. (51) is integrated to give 

x\r 1 )=x BB + 0{r l 2 ) . (252) 

The coefficients E a l [x 3 ' {rj)] are then Taylor-expanded: 

E a i [x^T,)]=E a i {x BB ) + 0(r 1 ) . (253) 

With expansions (118), (251) and (253), Eq. (51) is integrated to give an improved expansion 

x'W = *bb + K a Ej(xi B ) ry + 0( n 2 ) . (254) 
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F Fluid kinematical variables via boost transformations 

Let u denote a fluid 4-velocity vector field. Then fluid kinematical variables arise in the irreducible 
decomposition of the spacetime gradient of u, i.e., 

V a u b = -u a u b + Hh ab + a ab + u ab , (255) 

defining the fluid acceleration u a , the fluid Hubble scalar Ft, the fluid shear rate a ab , and the fluid vorticity 
Lu a b, respectively. 

With respect to an Eulcrian reference frame, comoving with a unit timclikc congruence u, we have 
u a =T(u a + v a ) , h a b = h a b + T 2 (v 2 u a + v a ) u b + T 2 (u a + v a ) v b , (256) 



where h a b :— S a b + u a u b , and the Lorentz factor is defined by T := 1/Vl — v 2 , with v 2 := v a v a . We 
define the peculiar kinematical variables 

9(v):=D a v a & ab (v):=-D {a v b) uj a (v) := e abc D b v c . (257) 

In the case when u is vorticity- free [i.e., uj a (u) = 0], we find the boost transformations: 

Fluid acceleration: 

h a b ii b = T 2 [ii a + v^ + (H+\6)v a + {a a b + a a b )v b + e abc £j b v c ] 

+ T 4 (v b v b + i 9v 2 + a bc v V) v a (258) 
u a u a = -u a v a . (259) 
Fluid Hubble scalar: 

H = T(H+l§+l u a v a ) + i T 3 (v a v a + \9v 2 + a ab v a v b ) . (260) 

Fluid shear rate: 

h c a h d b a cd = T(a ab + a ab )+T 3 v {a {u b) +v {b)) +a b)c v c ) + 2T 3 v (a 6- b)c v c 

- | T {ii c v c ) [ h ab + T 2 v a v b } - i T 3 {v c v c ) [ h ab - 2T 2 v a v b } 
+ i T 5 9 v a v b - i T 3 9v 2 [ h ab + T 2 v a v b ] - § r 3 (a cd v c v d ) [ h ab - 2T 2 v a v b ] (261) 
u a u b ~o ab = ~a ab v a v b (262) 
u b h c a o bc = - ~a bc v b h c a . (263) 
Fluid vorticity: 

h c a h d b Lv c d = Te abc Lo c + T 3 v [a (u b] +v {b)] + a b]c v c ) + 2T 3 v [a e b]cd u} c v d (264) 

u b h c a Co bc = -u bc v b h c a . (265) 

Equations (258), (260), (261) and (264) provide the generalizations to Go cosmologies of Eqs. (1.15), 
(1.16), (1.27) and (1.26) given by King and Ellis [18] for the SH subcase. 

The covariant derivatives of v that appear in Eqs. (258)-(265) convert into orthonormal frame ex- 
pressions according to 



v 



<«> 



e (v a ) -e af3 ~<rL v 7 (266) 

D a v a — (e a -2a a )(v a ) (267) 

D(a^fc) — > (e( Q + a< a ) (vp)) ~ eysianp) 1 v 5 (268) 

e abc D b v c — e"^ ( e/J - 0/ j) (v 7 ) - n a fi v fi . (269) 
Finally, we define Hubble-normalized fluid kinematical scalars by 

U 2 := U% )(*) = k (270) 



H J \H 

== = U^ ab ) (271) 



6 V H J \H 

UJab\ 

H )\H 



W 2 := U^](^-) = \{W ab W ab ). (272) 
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